High pressures make Hg a transition metal in a thermodynamically stable solid 
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The appropriateness of including Hg among the transition metals has been debated for a long 
time. Although the synthesis of HgF4 molecules in gas phase was reported before, the molecules 
show strong instabilities and dissociate. Therefore, the transition metal propensity of Hg remains 
an open question. Here, we propose that high pressure provides a controllable method for preparing 
unusual oxidation states of matter. Using an advanced structure search method based on first- 
principles electronic structure calculations, we predict that under high pressures, Hg can transfer 
the electrons in its outmost d shell to F atoms, thereby acting as a transition metal. Oxidation of 
Hg to the +4 state yielded thermodynamically stable molecular crystals consisting of HgF4 planar 
molecules, a typical geometry for d s metal centers. 

PACS numbers: 62.50.-p, 61.50.Ah, 61.50.Ks, 82.33.Pt, 81.40.Vw, 71.70.Fk 



A goal of physics and chemistry is to prepare unusual 
states of matter beyond the naturally occurring forms. 
One of the most challenging and attractive topics in this 
field is the preparation of high oxidation states of Hg. [1 , 
12] Group IIB elements, including Zn, Cd, and Hg are 
usually defined as post-transition metals [3] because they 
are commonly oxidized only to the +2 state. [U [5] Their d 
shells are completely filled and do not participate in the 
formation of chemical bonds; however, there is a high 
expectation that Hg should be stable at higher oxidation 
states due to the large relativistic effects that increase 
the energies of its 5c? levels. 

Earlier attempts focused on HgF 4 molecules in gas 
phases, including several quantum chemistry calcula- 
tions on the molecule. [6TTT3"] It was predicted that 
HgF 4 molecules can resist dissociation to HgF 2 and F 2 
molecules. [9] [14j [15] However, the thermodynamically 
stable form of HgF 2 at low temperatures is the solid 
phase rather than the gas phase. [2] HgF 4 molecules in 
the gas phase or in molecular crystals can, therefore, be- 
come unstable and dissociate to form HgF 2 and F 2 , as 
predicted by our calculations (see below). For this rea- 
son, the synthesis of HgF 4 molecules has shown to be ex- 
tremely difficult. Recently, after numerous experimental 
attempts, the trace of HgF 4 molecules was observed using 
matrix-isolation infrared (IR) spectroscopy. [1] The HgF 4 
molecules were unstable and dissociated to yield HgF 2 
and F 2 within very short periods of time. To demon- 
strate that Hg is truly a transition metal, one must ther- 
modynamically stabilize a Hg center in the +4 oxidation 
state. 

In this letter, we propose a totally new approach that 
utilizes high-pressure techniques to stabilize Hg 4+ in 
the solid phase. Over the past several decades, high- 
pressure methods based on diamond anvil cells (DACs) 



have yielded novel structures of elements and compounds 
under extremely high pressures (beyond megabars) . This 
field has been particularly advanced by the recent de- 
velopment of first-principles structure-searching methods 
that can predict high-pressure structures to extremely 
high accuracy without experimental or intuited input. 
For example, the recently implemented particle swarm 
optimization (PSO) technique [THJ E] successfully pre- 
dicted the high-pressure structures of various systems, 
including lithium [T5]. water ice[T2], and Bi 2 Te3.[20] In 
most high- pressure studies, the oxidation states of the 
compounds remain the same at high pressures. We show 
here that the high-pressure techniques have the potential 
to provide a controllable method for achieving unusual 
oxidation states. 

We conducted a structure search over the full structure 
configuration space, based on ab initio total-energy cal- 
culations and a PSO technique, implemented in the CA- 
LYPSO (crystal structure analysis by particle swarm op- 
timization) code. [161 117] The underlying ab initio struc- 
tural relaxations and electronic band structure calcu- 
lations were performed using the framework of den- 
sity functional theory (DFT), implemented in the VASP 
code.[H] The generalized gradient approximation (GGA) 
in the framework of Perdew-Burke-Ernzerhof (PBE) [22 
was used to treat the exchange-correlation functional, 
and the projector augmented wave (PAW) [23] pseudopo- 
tentials were used to describe the ionic potentials. The 
cutoff energy for the expansion of the wave function into 
plane waves was set to 700 eV, and the Monkhorst-Pack 
k meshes [53] were chosen to ensure that all enthalpy 
calculations converged to better than 1 meV/atom. The 
structure searches were performed with a unit cell con- 
taining 20 atoms (4 Hg and 16 F) and were performed 
over the pressure range 0-300 GPa. The phonon spectra 
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FIG. 1: (Color figures available online.) (a) and (b): Two views of 
the Pca2\ structure, the optimized lattice parameters at 15 GPa 
are a = 14.183 , b = 3.478 and c = 4.153 with Hg occupying 
4a (0.082, 0.277, 0.486) and F occupying 4a (0.309, 0.530, 0.321), 
4a (0.306, 0.173, 0.523), 4a (0.663, 0.211, 0.337) and 4a (0.521, 
0.771, 0.982) positions; (c) side view and (d) top view of the I4/m 
structure, the optimized lattice parameters at 150 GPa are a = 
4.366 and c = 3.498 with Hg occupying 2b (0, 0, 0.500) and F 
occupying 8h (0.304, 0.888, 0) positions; (e) top view and (f) side 
view of the IA/mmm structure; (g) the CaF2 structure; (h) the 
HgCb structure; (i) the F2 molecular crystal structure. The large 
dark-blue balls and the small light-grey balls indicate Hg and F 
atoms, respectively. Parameters of more related structures can be 
found in the online supplementary materials. 
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FIG. 2: (Color figures available online.) (a) Comparison of the 
enthalpies as a function of the pressure for the Hgi-F4 structures 
selected from the PSO search. The results are reported with re- 
spect to the enthalpy of the Pca2\ structure, which is the most 
stable structure under ambient conditions. (b) Energy of the 
HgF2(s)+F2(s)— ^ HgF4 (s) reaction. The enthalpy of H5F4 was 
obtained from the 74/m structure. The blue dashed line shows the 
enthalpy gain due to the volume reduction, PAV . The area shaded 
in light blue indicates the range of pressures under which HgF4 is 
thermodynamically stable at K. 



were calculated using the Phonopy program431 

As shown in Fig. [2ja), two thermodynamically sta- 
ble structures were identified over the pressure range of 
0-300 GPa. The most stable structure under ambient 
conditions had a symmetry of Pca2\ and consisted of 
chain- like HgF2 and F2 molecules [FigjTJa) and (b)], in- 
dicating that at low pressures, the Hgi-F4 system tended 
to decompose into HgF2 and F2. The shortest F-F dis- 
tances were 1.485 A at GPa, similar to the F-F bond 
length in the gas phase. [35] At a pressure of 20 GPa, the 
system transformed into a structure with IA/m symme- 
try [Figjljc) and (d)]. This structure consisted of HgF4 
molecules and each Hg had four nearest neighboring F 
atoms. At 50 GPa, the Hg-F distance was 1.949 A. This 
distance was significantly shorter than the average Hg-F 
distance of 2.215 A at the same pressure. On the other 
hand, the shortest F-F distance was 2.112 A which was 
significantly longer than the F-F molecular bond length 
of 1.439 A, indicating that the molecular features of F2 
were destroyed upon formation of the HgF4 molecules. 
The structural parameters for both structures and the 
dependence of bond lengths in IA/m structure on in- 
creasing pressure can be found in the online supplemental 
material. 

Another structure displayed a symmetry of I4/mmm 



[Fig. ITje) and[TJf)]. Although this structure appeared 
similar to the 74/m structure as viewed from the top, 
each Hg atom was bonded to eight F atoms to form a 
chain. The enthalpy of this structure was lower than that 
of the Pca2\ structure at 100 GPa; however, over the 
pressure range 0-300 GPa, the enthalpy exceeded that of 
the IA/m structure, a HgF 4 molecular crystal. 

As shown experimentally[5Sl [57] and in our calcula- 
tions, HgF2 was stable in extended solid structures other 
than the HgF2 chains present in the Pca2\ structure of 
the Hgi-F4 system. To determine whether the separated 
HgF 2 and F 2 phases could react to form a HgF4 phase, we 
calculated the energy of the HgF 2 (s)+F 2 (s)— !• HgF 4 (s) 
reaction. The structures of HgF 2 and F 2 under pressure 
were searched using the PSO method. HgF2 was found to 
be stable in the CaF 2 structure at to 100 GPa,[26l|27] in 
the low-pressure structure of HgCl 2 (Prima) [2"5] from 100 
to 300 GPa. We also found that F2 maintained its molec- 
ular features and was stable in its low-pressure structure 
C12/clg2] up to 300 GPa. 

As shown in Fig. [2](b) , our results clearly showed that 
at ambient pressures, the separated HgF 2 and F 2 phases 
were stable against the formation of HgF 4 , i.e. the reac- 
tion energy was positive. The reaction energy decreased 
rapidly as the pressure increased and becomes negative 
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FIG. 3: (Color figures available online.) (a) Phonon spectrum of 
the HgF4 molecular crystal in the IA/m structure under a pres- 
sure of 150 GPa. The energy regions that corresponded to the 
F in-plane, the F out-of-plane, and the Hg motions are shaded 
in red, green, and white, (b) The phonon projected density of 
states (PDOS) of HgF4 in the IA/m structure at a pressure of 150 
GPa. The black line and the red and green shaded areas show the 
DOS projected onto the Hg, the F in-plane, and the F out-of-plane 
motions, (c) The summed phonon PDOS of HgF2 in the Prima 
structure and F2 in the C12/cl structure at 150 GPa. The black 
line shows the motions of Hg in HgF2. The blue and purple shaded 
areas show the motions of F in the HgF2 and F2 solids, respec- 
tively, (d) The three selected IA/m phonon frequencies at T as a 
function of the pressure. The black, orange, and blue lines denote 
the Hg-F stretching and bending modes, as shown in the insets, (e) 
The calculated reaction free-energies as a function of the pressure 
at 0, 1000, 2000, 3000, and 4000 Kelvin. 



at 60 GPa, indicating that HgF4 became thermodynami- 
cally stable. The energy gain upon formation of the HgF4 
molecules was quite large and increased significantly with 
the increasing pressure. At a pressure of 150 GPa, the 
energy gain was 0.8 eV or 77 kJ/mole, providing a large 
driving force for oxidation of Hg to the +4 state. 

The phonon spectrum of HgF4 in the IA/m structure, 
of HgF 2 in the Pnma structure, and of F 2 in the C12/cl 
structure were calculated over the pressure range of 0- 
300 GPa. The results show that the structures of all three 
solids were dynamically stable up to 300 GPa, i.e. none 
of the phonon modes featured imaginary frequencies. As 
shown in Fig. [3]ja) and (b), the low-frequency modes of 
the HgF4 molecular crystal at 150 GPa (from to 10 
THz) were dominated by the motions of the Hg atoms, 
whereas the modes in the frequency range 10-18 THz 
and 18-25 THz were dominated by the F out-of-plane 
(z) and in-plane (x and y) motions. The Hg-F stretch- 
ing modes in HgF4 reached their peaks at 24 THz (800 
cm -1 ). On the other hand, the Hg-F stretching modes in 
HgF 2 reached their peak around 20 THz (650 cm -1 ), as 
shown in Fig. [3^c) . The F-F stretching modes were much 



higher in energy, and a frequency gap was observed from 
21 to 31 THz (700-1000 cm" 1 ) in a HgF 2 and F 2 mixed 
sample. The vibrational peak observed in this gap will 
indicate the formation of HgF4 molecular crystals. The 
HgF4 vibrational modes in this frequency gap included 
Hg-F stretching and bending modes and could be either 
IR or Raman active. As shown in Fig. [3](d), the frequen- 
cies of these peaks tended to split under higher pressures, 
indicating a broadening of the HgF 4 signature peaks with 
increasing pressure. 

Another important issue associated with the synthesis 
of HgF4 in a high-pressure DAC is the persistence of the 
stable structure at high temperatures. Although HgF 4 is 
thermodynamically stable at K, it does not form easily 
at very low temperatures in a DAC because of the large 
kinetic barrier to the solid-state reaction. To overcome 
this barrier, thermo-heating or laser-heating techniques 
may be used to increase the temperature in the DAC 
up to a few thousand degrees Kelvin. We calculated the 
reaction energies at finite temperatures using the calcu- 
lated phonon spectra of the HgF4, HgF 2 , and F 2 solids. 
As shown in Fig. |3je), elevated temperatures tended to 
destabilize the HgF4 solid; however, even at T = 4000 K, 
HgF 4 could still reach thermodynamic stability at pres- 
sures exceeding 75 GPa. 
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FIG. 4: (Color figures available online.) (a) The crystal field split- 
ting diagram of a square planar complex; the light-blue balls in- 
dicate the F anions and the orange and red spheres indicate the 
d x 2_ y 2 orbit; (b) the d orbitals were split under the square pla- 
nar crystal field; (c) the calculated electron localization function 
(ELF) for HgF 4 in the IA/m structure at 150 GPa; (d) the calcu- 
lated PDOS for HgF 2 in the HgCb structure at 50 GPa and for the 
HgF4 molecular crystal in the IA/m structure at 50 and 150 GPa. 
In the PDOS plots, the black solid lines indicate the total DOS, the 
blue solid lines indicate the F 2p states, the red and green shaded 
areas indicate the Hg 6s and 5d states. 
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The planar geometry of the HgF4 molecules is typi- 
cal of d 8 configurations, indicating the involvement of 5c? 
electrons in the formation of Hg-F bonds. As shown in 
Fig. EJa) and (b), the d orbitals split into four groups 
in a square planar crystal field due to differences in the 
interactions with the four surrounding F anions. Among 
these orbitals, d x 2_ y 2 was the highest in energy. The 
transition metal cations could be energetically stabilized 
in this configuration in the event that 8 electrons filled 
the 4 lower d states and left the d x 2_ y 2 states empty. The 
bonding features of HgF4 were also examined by calculat- 
ing the electron localization function (ELF) [35] using the 
VASP code. The results [Fig. [2fc)] clearly showed the 
equivalent covalent bonds of Hg and its four neighboring 
F atoms. 

More direct evidence for the involvement of d elec- 
trons in the formation of chemical bonds with F was 
obtained from the projected density of states (PDOS). 
We calculated and compared the total DOS and PDOS 
for HgF4 in the 14/ m structure and HgF2 in the HgCi2 
structure at 50 GPa. HgF 2 revealed a gap of 1.5 eV. The 
conduction bands consisted mainly of the F p and Hg 6s 
states, whereas the valences bands consisted of both F p 
and Hg 5o! states [Fig. gd)]. The PDOS showed that 
all 5d states were filled, and Hg was apparently in the 
+2 state. On the other hand, HgF4 displayed a much 
smaller gap of only 0.71 eV. More importantly, the 5d 
states were distributed over both the valence and the 
conduction bands, indicating that the 5d electrons were 
depleted, in agreement with the planar geometry of the 
HgF4 molecules and the d s configuration analysis. 

Several effects of the high pressure assist in stabiliz- 
ing the HgF4 molecules. One mechanism is based on 
the crystal field splitting. Because high pressures re- 
duce the inter-atomic distances, including those between 
Hg and its neighboring F atoms, the high pressures en- 
hance the crystal field splitting. As a result, the energy 
of the d x 2_ y 2 state increases with increasing pressure. 
The PDOS of HgF 4 in the 14/ m structure at 150 GPa is 
shown in Fig. |4jd). A comparison with the PDOS at 50 
GPa shows that although the gap changed only slightly 
with increasing pressure, the center of the F 2p states 
was significantly lower relative to the d x 2_ y 2 state, indi- 
cating a larger energy gain with electron transfer from 
the bd to the 2p orbital at higher pressures. Another im- 
portant factor is the reduction of the volume during the 
formation of solid HgF4. This factor further increases the 
enthalpy gain. As shown in Fig. gb), the value of PAV 
increased with increasing pressure. Although its contri- 
bution to the enthalpy gain was small at low pressures, 
this contribution became dominant at pressures exceed- 
ing 50 GPa. The reduction in the volume arose from the 
shortening of Hg-F bonds as Hg changed from the +2 to 
the +4 oxidation states due to the involvement of the 5d 
electrons in the formation of stronger Hg-F bonds with a 



larger covalent component. When the enthalpy gain was 
sufficiently large, the reaction HgF 2 (s)+F 2 (s) —> HgF4 
(s) became exothermic and HgF4 became thermodynam- 
ically stable. As calculated, this occurred at 60 GPa. 

In summary, we demonstrated that under an external 
pressure applied using the DAC technique, Hg could 
be thermodynamically stabilized in a high oxidation 
state to form solid structures of HgF4 molecules. The 
atomic structure and the electronic structure of the high- 
pressure phases suggested that the 5d orbitals of Hg were 
involved in chemical bonding. Thus, our calculations 
suggested that Hg behaves as a true transition metal 
under high pressures. High-pressure techniques have 
been used to alter the chemical inertness of noble gases 
and to sustain them in unusual chemical states. [55H55] 
Our study also suggests that high-pressure techniques 
may be used as a controllable method to achieve unusual 
oxidation states in matter. 
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